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Abstract 



r*| ■ The theory of matchgates is of interest in various areas in physics and 

C^ ' computer science. Matchgates occur in e.g. the study of fermions and 

? , , spin chains, in the theory of holographic algorithms and in several recent 

'^Z works in quantum computation. In this paper we completely characterize 

the class of boolean functions computable by unitary two-qubit match- 
es] ■ gate circuits with some probability of success. We show that this class 
^ ' precisely coincides with that of the linear threshold gates. The latter is a 
CO I fundamental family which appears in several fields, such as the study of 
^T ■ neural networks. Using the above characterization, we further show that 

the power of matchgate circuits is surprisingly trivial in those cases where 
the computation is to succeed with high probability. In particular, the 
1/^ ■ only functions that are matchgate-computable with success probability 

i~} ' greater than | are functions depending on only a single bit of the input. 

O' 

1 Introduction 

^\f ' One of the great virtues of the field of quantum computation is that it intercon- 

3 , nects fundamental questions in physics and computer science. The concept of 

the quantum computer [1] precisely captures the intrinsic computational power 
locked within quantum mechanics |2], and makes it possible to address deep 
problems such as the relationship between quantum and classical computational 
capabilities [IIIHIS]. At the same time, within the theory of quantum computa- 
tion it is possible to characterize, in a precise sense, how "hard" is it to simulate 
physical systems of interest, such as certain ground state problems [6j and time 
evolutions [7,. 

Of particular interest, also in recent work, is the class of quantum processes 
generated by matchgates O [3 [9l [TOl HH [l2l US] . The latter are a class of uni- 
tary two-qubit operations that are defined by certain algebraic constraints. The 



theory of matchgates is an instance of a research area that displays strong con- 
nections to both physics and computer science [5j [8j [9l [101 HH Ull US] . In the 
study of strongly correlated systems, for example, the dynamics of an important 
class of ID quantum systems such as the XY model are modeled by matchgate 
circuits i.e. for such hamiltonians H one can construct a poly-size matchgate 
circuit Ct such that Ct = e**^ for any time t (see e.g. [TO])- Employing mappings 
between spin-^ systems and fermions, matchgate circuits further describe the 
dynamics of all non- interacting fermionic systems [5] . In the theory of quantum 
computation, matchgates are of particular interest as they provide a key exam- 
ple of class of nontrivial quantum circuits that cannot offer any speed-up over 
classical computers (in spite of e.g. the complex entangled states such circuits 
may generate) [10] . In addition, matchgate computations were recently found 
to be equivalent to space-bounded universal quantum computation [12] . In clas- 
sical computer science, finally, matchgates occur in various studies related to 
e.g. the theory of holographic algorithms [SKH]. 

The aim of the present paper is to characterize the computational power 
of matchgate circuits. We will in particular study which boolean functionqj 
can be computed with such circuits. Given an arbitrary matchgate circuit U, 
the question is asked which boolean function f{x) can be computed (proba- 
bilistically) by initializing the system in the computational basis state |a;,0) 
(where x represents an input string and is a string of ancillary zeros), by 
subsequently running the circuit U and finally measuring, say, the first qubit. 
This setting captures in perhaps the most elementary way the computational 
power of matchgate circuits, associating which each circuit a yes/no question as 
commonly done. 

We remark that, beyond its natural computer scientific interest, such an in- 
vestigation is relevant from a intrinsic physical perspective as well. In particular, 
given the aforementioned equivalences between matchgate circuits, fermionic 
systems and ID spin systems, the present work aims at gaining insight in the 
link between the physics of these systems and their computational capabilities. 
In this context one may pose a variety of interesting questions such as: 'Does the 
presence of a quantum phase transition in the XY model leave any signature in 
the associated class of functions which can be computed by time-evolving such 
systems?' The present work is also situated within such a program. 

In the following we will characterize the family of matchgate-computable 
functions in full generality. We will find that this class precisely coincides with 
the class of linear threshold gates (LTGs) |14] . The latter is a fundamental 
family of functions that has been a topic of study since the 1960s and that 
plays an important role in numerous areas. LTGs occur in the study of neural 
networks where these functions serve as elementary models of neurons |15| and 
in circuit complexity theory |16| : cf. also e.g. |17| and references within for a 
number of recent investigations on LTGs. The existence of a connection between 
matchgates and LTGs may be considered surprising, since a priori there is no 



^Henceforth, whenever we refer to a 'function' we will often mean a boolean function. This 
will be clear from the context. 



obvious relation between these two theories. 

Below we state the contributions of this work in more precise terms. Here 
we highlight some particularly noteworthy aspects of our results. 

An interesting phenomenon occurs when considering those functions that are 
matchgate-computable with high success probability. Since a generic matchgate 
circuit may be a rather complex object, one would expect that this class of func- 
tions has some nontrivial character as well. We will show that this intuition is 
incorrect: any function that is matchgate-computable with success probability 
greater than 3/4 is proved to be trivial in that any such function can only de- 
pend on a single bit of the input. The origin of this apparent paradox is the 
strong set of constraints that are placed on any circuit (matchgate or conven- 
tional) that is to compute the correct function value with high probability on 
all inputs. Indeed, a generic circuit does typically not fulfill these requirements 
and hence does not meaningfully compute any function. The present example 
thus highlights the significance of bounded-error constraints in a rather striking 
way: in spite of the potentially complex structure of general matchgate circuits, 
the only instances which turn out to satisfy the bounded-error constraints are 
circuits computing mere single-bit — i.e. utterly trivial — functions! 

Note that the above feature is of interest from a physics perspective as well: 
it shows that the nontrivial physical properties (such as e.g. the presence of a 
quantum phase transition) of a family of quantum systems are not guaranteed to 
translate into any nontrivial associated computational model. Indeed, in spite of 
the interesting physical processes modeled by matchgates, these processes turn 
out to have trivial power when used as a computer which is to solve problems 
with high success probability. 

Overall, the results in this paper indicate that matchgate circuits have a 
rather weak computational power. This is in part reflected in the above phe- 
nomenon, but also in an additional result obtained in this work. We will char- 
acterizes the power of matchgate circuits in terms of a strikingly elementary 
classical computer which is capable of computing every matchgate-computable 
function with the same success probability as the optimal matchgate circuit. 

Notation. — In the following [n] denotes the set of positive integers from 

1 to n, for any n. If x = (xi, . . . ,a;„) is a string of bits Xk S {0, 1}, then x 
denotes the ±1 vector that is obtained by replacing all 0-components of the n- 
bit string x by 1 and all 1-components by —1. The symbol ® denotes addition 
modulo 2. The 1-norm of a complex vector w = (wi, . . . , Wn) is denoted by 

Iklli — Efckfel- 

2 Main results 

Here we summarize the contributions of this paper. To do so, we state some 
preliminary definitions. 



A matchgatcl G is any two-qubit operation with matrix representation 
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in the standard basis, where A and B belong to SU{2). We will consider cir- 
cuits composed of matchgates acting on nearest-neighbor qubits (assuming a 
one-dimensional ordering of the qubits) i.e. the standard scenario in which 
matchgates are considered. 

Let / : {0, 1}" —i' {0, 1} be a boolean function on n bits and let U be an m- 
qubit unitary operation with m > n. We say that U computes / with probability 
at least p if for every n-bit string x, the preparation of the m-qubit state U\x, 0) 
(where denotes a string of m — n zeroes) followed by a computational basis 
measurement of the first qubit yields the outcome f{x) with probability at 
least p. An n-bit boolean function / is said to be matchgate-computable with 
probability at least p if there exists a matchgate circuit on m qubits, for some 
m > n, which computes / with probability at least p. Note that in the latter 
definition no restriction is placed on the size of m as compared to n, nor on 
the number of gates in the matchgate circuit compared to n. However, below 
we will find that every matchgate-computable function can be computed by 
a matchgate circuit acting on at most n -I- 1 qubits (without decreasing the 
success probability). Moreover, it is known that any matchgate circuit acting 
on m qubits can be re-rexpressed as a matchgate circuit of size 0{m^) (see |10j). 



I. Main Theorem. — We will show that the class of matchgate-computable 
functions coincides with the family of linear threshold gates. A boolean function 
f on n bits is called a linear threshold gate (LTG) if there exist and an n- 
dimensional real vector w and a real constant 6 such that f{x) equals if and 
only if w'^x + 9 is strictly positive. The vector {w, 9) is called a representation 
of /. Examples of linear threshold gates are the NOT gate, the n-bit AND and 
OR and the majority gate. We will consider linear threshold gates that are 
supplemented with a parameter which is called the margin of the LTG. Given 
a linear threshold gate / on n bits with representation {w,9), the margin of 
this representation is defined to be the minimal value of |w"^x-|-0| over all n-bit 
strings x. The margin e of / itself is the maximal achievable margin of any 
representation {w, 9) of / which is normalized in the sense that ||w||i -I- \9\ = 1. 

The main result of this paper achieves a complete characterization of all 
matchgate-computable functions: 

Theorem 1. Let f be a boolean function on n hits and let p G (0.5, 1]. Then 
the following statements are equivalent: 



^The term 'matchgate' sometimes refers to a larger class of operations (tensors) which 
may be both non-unitary and which may act on more than two qubits, containing the unitary 
two-qubit gates ((!} as a subclass; cf. lOl llll . In this paper, however, a matchgate is always 
taken to be a unitary two-qubit operation as defined in ifT) . 



(a) f is matchgate- computable with probability at least p. 

(b) f is a linear threshold gate with margin e > 2p — 1. 

Moreover, (a) holds if and only if there exists a matchgate circuit acting on at 
most n + 1 qubits which computes f with probability at least p. 

Note that theorem [1] connects the margin e of an LTG with the optimal success 
probabihty p of computing this function using matchgate circuits. In particular, 
e is small iff p is small. This implies that the class of linear threshold gates 
computable with matchgate circuits grows larger as the required probability of 
success is decreased. When p is allowed to be arbitrary close to 0.5, the full 
class of LTGs is matchgate-computable due to theorem [TJ 

As an immediate corollary of the above result, it follows that matchgate 
circuits do not have universal classical computational power, even when an 
unbounded error is allowed i.e. a success probability strictly greater then 0.5 
but (with increasing n) potentially exponentially close to 0.5. This property 
follows immediately from the elementary fact that there exist functions that are 
not LTGs, such as the two-bit parity gate. 

II. Bounded-error computations. — Surprisingly, it follows from theo- 
rem [T] that matchgate circuits can only compute trivial functions if the compu- 
tation is to succeed with high probability: 

Theorem 2. A boolean function f is matchgate-computable with probability 
p > J if and only if this function is either constant or depends on a single bit 
of its input. 

Due to theorem m any function that is matchgate-computable with probability 
p > 3/4 is a linear threshold gate with margin e > 1/2. We will show that 
the only LTGs having such large margin are constant or depend on one input 
bit, leading to the proof of theorem [21 Note that functions which depend on 
a single input bit have the very simple form f{x) = Xk or f{x) — 1 — Xk for 
some fc, i.e. / returns the fc-th bit of its input or its negation. As discussed 
in the introduction, theorem [5] is a somewhat unexpected result, given that 
generic matchgate circuits are rather nontrivial objects which e.g. describe 
physical systems that may exhibit interesting behavior such as e.g. quantum 
phase transitions. In spite of this rich structure, theorem [2] shows that the 
computational power of matchgate circuits is near- vanishing in those cases where 
the correct answer is to be produced with high probability. 

Of particular interest are computations which yield the correct output with 
a probability that is bounded away from 0.5 by an inverse polynomial in the 
input size. Let J- — {/„ : n = 1, 2, . . . } be a family of boolean functions where 
/„ acts on n bits and let {p„} be a family of probabilities where p„ > 0.5 and, 
for n large, Pn is bounded away from 0.5 by an inverse polynomial in n. We 
say that J-" is matchgate-computable with poly-bounded error if there exists 
such a poly-bounded family of probabilities as well as a family of matchgate 
circuits {Un}, such that C/„ computes /„ with probability at least p„. It is 



standard that any computation with poly-bounded error can be promoted to an 
almost-deterministic computatiorlf] by repeating the computation poly(n) times 
and taking the majority vote of all obtained results. Due to theorem [U the 
families of boolean functions that are matchgate-computable with poly-bounded 
error precisely coincides with those LTGs having poly-bounded margin i.e. the 
margin e„ of /„ scales an an inverse polynomial in n. We will furthermore 
show that a family of LTGs has poly-bounded margin if and only if /„ has a 
representation (wmOn) where the coefficients of Wn and 6n are integers that are 
at most polynomially large; such families of threshold gates are said to have 
polynomial integer weight. This leads to the following concise characterization. 

Theorem 3. A family of boolean functions is matchgate-computable with poly- 
bounded error if and only if it is a family of linear threshold gates with polynomial 
integer weight. 

Remark. — The near-determinstic computation of n-bit LTGs with polyno- 
mial integer weight is obtained by running the computation poly(n) times and 
computing the majority vote of all measurement outcomes. It is intriguing that 
the majority function is indeed an LTG — and hence matchgate-computable — 
however in order to properly amplify the success probability, the majority gate 
itself needs to be computed with suitably high success probability (the latter 
e.g. being exponentially close to 1). However, due to theorem [2] the majority 
function cannot be computed with good success probability by any matchgate 
circuit. Therefore, to obtain the proper amplification a final non-matchgate 
computation is needed to compute the majority vote — even though the latter 
comes intriguingly close to being suitably matchgate-computable! o 

III. An equivalent classical computer. — As a final result, we will con- 
struct a classical computational scheme that is equivalent to the matchgate cir- 
cuit model in the following sense: any function which is matchgate-computable 
with probability at least p can be computed with probability at least p with 
our classical scheme, and vice versa. As we will show, the required classical 
computer is very simple, as it essentially requires a single sample of a fixed (i.e. 
independent of the input) probability distribution on the set of integers from 1 
to n-|- 1, together with the possibility of performing a single bit flip of one of the 
input bits, depending on the outcome of the sampling. This characterization is 
a further illustration (in addition to e.g. theorem [5]) of the weak computational 
capabilities of matchgate circuits. 

Whereas the general result will be stated below, here we illustrate the scheme 
with an example. Let n be odd and consider the majority function /maj on n- 
bits, which is an LTG. It can be shown that the margin of fmaj is Cmaj = n~^. 
Due to theorem [I] the optimal success probability of computing /maj with a 
matchgate circuit is 

Pmaj = -^ = 2 + 2^^- ^^^ 



^That is, the success probability is exponentially (in n) close to 1. 



Now consider the following elementary classical computation: 

• Choose an n-bit input string x; 

• Generate a random integer k between 1 and n; 

• Output the bit value Xk- 

It can easily be shown that, for every x, the output of this computation is 
/maj(a^) with probability at least Pmaj- In other words, the above elementary 
classical computation is capable of computing the majority function with the 
same probability of success as the optimal matchgate circuit can! 

The above example is not coincidental, as we will show that every LTG / 
can similarly be associated with a simple classical computation of the above 
nature. 



3 Matchgates 

In this section we recall some basic properties of matchgates, which were defined 
in Eq. ([T]). We emphasize that, henceforth, the term 'matchgate circuit' will 
always refer to a circuit composed of matchgates acting on nearest-neighbor 
qubit lines, as commonly done. 

We first recall a celebrated result about the classical simulation of matchgate 
circuits, first proved in |S] and subsequently investigated by a series of other 
authors [i[i[ini[n]. 

Theorem 4. Consider a poly-size n-qubit circuit composed of matchgates acting 
on nearest-neighbor qubits. The circuit acts on an arbitrary standard basis input 
and is followed by a standard basis measurement of the first qubit. Let po denote 
the probability of obtaining the bit as an outcome. Then there exists a classical 
algorithm which computes po up to m bits in poly{n, m) time. 

As a consequence of this result, there exists a poly-time classical algorithm 
to sample from (a distribution that is exponentially close to) the distribution 
{po, 1 — Po} i-G. any matchgate computation of the above type can be simulated 
classically in poly-time. 

We discuss some further well-established properties of matchgates, which 
will be used in the proof of theorem [T] We refer to e.g. [TD] for elementary 
proofs of these properties. 

First, without loss of generality, in the following we will always consider poly- 
size families of matchgate circuits, as it is known that any, possibly exponential 
size, n-qubit matchgate circuit family can be re-expressed a matchgate circuit 
family of size 0{n^). Second, consider the n-qubit Jordan-wigner operators: 

fc-i fc-i 

C2fc-i = -'^fe J]^ Zi, C2k = ^fc n ^'' ^"^) 



where X^, Yfc and Zk denote the Pauh X, Y and Z operators acting on qubit 
fc, respectively. Then an n-qubit unitary operation C/ is a matchgate circuit (up 
to a global phase) if and only if there exists an operator R G SO{2n) such that, 
for every ^ £ [2n], it is the case that 



2n 

C/tc^[/ = ^i?^,c,. (4) 

It further holds that [/ is a matchgate circuit if and only if there exists a her- 
mitian operator H lying in the linear span of the products c^Ci, (where ^ ^ v), 
such that U oc e*^ . Such an operator H is sometimes called a quadratic hamilto- 
nian. Quadratic hamiltonians describe the physics of systems of non-interacting 
fermions. We will not discuss this connection to fermionic physics here as it 
would lead us too far outside of the scope of this work, and we refer to e.g. [5]. 
The following are some examples of matchgates and matchgate circuits. The 
fermionic SWAP (fSWAP) operation is easily seen to be a matchgate. This 
operation sends the basis state \ab) to {—l)°'^\ba), for every a,b = 0,1, i.e. it 
swaps the qubits and adds an overall minus sign if both qubits are in the state 1 1} . 
Other elementary examples of matchgate circuits are the products Cfj_Ci, for any 
jjL ^ V. Indeed, any such product can be written as an exponential of a quadratic 
hamiltonian. Denoting H :— iCfj,Ci, it is easily verified that H = H^ = H^. 
Hence 

e"^" (X H (X Cfj^c^, (5) 

where in the first identity we have used that e**'^ = (cosi)/ + i{smt)H since 
H = H^. To obtain a more nontrivial example of a matchgate circuit, consider 
the hamiltonian of the one-dimensional transverse ising model: 

n — 1 n 

Hi^-Y, JkXkXk+1 - J2 hkZk, (6) 

fe=l k=l 

for some real constants Jj, and h^. It can readily be verified that 

XkXk+i (X C2kC2{k+i)-i and Zk (x C2k-iC2k, (7) 

showing that Hi is a quadratic Hamiltonian. Therefore, for every real t the 
time evolution operator e''^^ can be written as a (poly-size) matchgate circuit. 
The notion of a matchgate-computable boolean function was introduced in 
section[2] Here we discuss some simple examples. It can easily be seen that every 
boolean function which is either constant or which depends on a single input 
bit can computed with unit probability by an elementary matchgate circuit 
composed of operators c^j,Ci, and fSWAP gates, acting on n -I- 1 qubits. Note 
that the only possible functions of this type are the functions a; — >■ 0, a: —?> 1, 
X ^ Xk and a; — s> \ — Xk, for some k. We show that a; -» 1 — Xfc is matchgate- 
computable with unit probability; the other three cases are treated similarly. 



For each k € [n] , denote the operator Uk acting on n + 1 qubits by 

n 
Uk '■— C2/c-lC2(ri+l)-l OC XkXn+1 | J_ Zj. (8) 

Note that Uk\x,0) oc |a;^fc,l), where x-,fc denotes the bit string obtained by 
flipping the fc-th bit of x. Now consider an elementary matchgate computation 
where first the state Uk\x,0) is prepared, followed by measurement of qubit 
k. This computation yields the bit 1 — Xk with unit probability. This shows 
that the desired function can be computed with unit probability by applying 
a suitable matchgate circuit and measuring some qubit in the computational 
basis. By applying a suitable sequence of fSWAP gates immediately before the 
measurement, we may assume w.l.o.g. that the first qubit is measured. Indeed, 
as the measurement is in the computational basis, the minus sign of the fSWAP 
gate has no relevance, and this gate acts as a simple SWAP. 

Remarkably, as we will prove in theorem [5J constant and single-bit func- 
tions are the only functions that are matchgate-computable with high success 
probability. 

4 Linear threshold gates 

In this section we discuss some elementary features of linear threshold gates 
(LTGs). For convenience we recall here their definition. An n-hit boolean 
function / is a linear threshold gate if there exists an n-dimensional real vector 
w and a real number 9 such that 

(-l)^(^) ==sign(w^i:H-6l) (9) 

for every n-bit string x. The pair {w, 6) is called a representation of /. Definition 
^ has an elementary geometrical interpretation. Taking an arbitrary w and 
9, the linear equation w^ z + = defines a hyperplane which divides the n- 
dimensional real space in two parts, say H^ and 7J_, where H^ consists of all 
z e M" such that vo^ z -\- > Q and _ff_ is defined as the complement of i/+. 
The LTG / associated with [w, 9) then simply evaluates whether a given {±1}- 
vector X lies in H^ or H^. Stated differently, a boolean function / is a linear 
threshold gate iff there exists a hyperplane in n-dimensional real space which 
separates the sets of inputs x that are mapped to and 1, respectively. 

It can easily be verified that the constant and single-bit functions are linear 
threshold gates. Other examples are the the n-bit AND, OR and majority 
function. The AND gate, for example, has a representation {w, 9) given by 
w — (!,...,!) and 9 = —n + ^. It is also known that not all functions are 
LTGs. For example, the two-bit parity gate 

f{xi,X2) ■■= Xi ®X2 (10) 

is not a linear threshold gate, as can be easily verified. 



Every linear threshold gate has infinitely many representations. For example, 
rescaling the vector (u;, 9) with a positive multiplicative constant trivially leads 
to the same associated function. We will say that the representation {w,6) 
is normalized if the 1-norm of this vector is equal to 1, i.e. ||w||i + \0\ = 1. 
Interestingly, every linear threshold gate has a representation {w, 9) where each 
Wk and 9 are integers. The intuition behind this result is the following. It is 
easy to show that every LTG has a representation (w, 9) where each Wk and 9 
are rational numbers, say Wfe = f^ and 9 = ^ for suitable integers ak,hk,c and 
d] this essentially follows from the fact that the rationals are dense in the reals 
and the property that small perturbations of a representation do not change the 
associated LTG. Multiplying the rational representation {w, 9) with the product 
in l^fcllMI then yields an integer representation. 

Further, it is known that every linear threshold gate on n bits has a repre- 
sentation (u", 9) where each of the components of w, as well as the number 0, are 
integers not greater than 2''^ in absolute value, with N = 0(n log n) [IH]; more- 
over there exist linear threshold gates where such large numbers are required 
[T5] . This shows that every linear threshold gate / admits a representation 
that can be fully described in terms of 0{ii? logn) bits, and that f[x) can be 
evaluated in poly-time when this particular representation is provided. 

Finally, we recall that that, given a representation {w,9) of an LTG /, the 
margin of this representation is defined to be the minimal value of \viFx + 9\ 
over all n-bit strings x. The margin e(/) = e of / is then the maximal margin 
over all normalized representations. Note that e(/), which should in principle 
be defined as a supremum, is indeed a maximum. This can be argued with 
standard methodo In section [6] we will focus in more detail on the properties 
of margins of linear threshold gates, which will lead to the proofs of theorems [2] 
and El Before doing so, we provide the proof of theorem [1] in the next section. 

5 Proof of theorem [1] 

In this section we prove theorem[T] The proof will proceed in three steps. In step 
1 we reduce the study of matchgate-computable functions to the investigation 
of matrix elements of the form {x\U^ ZiU\x) , where \x) denotes a computational 
basis state, f7 is a matchgate circuit and Zi is the Pauli Z operator acting on 
the first qubit. In step 2, which represents the main ingredient of the proof 
of theorem (TJ the most general form of such matrix elements is characterized. 
Finally, in step 3 the proof is completed by combining steps 1 and 2. 



*To see this, first note that the margin of a representation {w,9) = 2 is a, continuous 
function of z, being defined as the minimum of a finite number (i.e. 2" ) of continuous functions 
ga;(z) := \w'^x + d\ for every x. Thus e(/) is the supremum of a continuous function in z, over 
all representations z of / with ||z|ji = 1. As the set of all such normalized representations is 
a compact set, it follows that the supremum is reached. 
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5.1 Step 1 

Let / be an n-bit boolean function, let U be an arbitrary ?Tj-qubit unitary 
operation for some m> n and fix p e (0.5, 1]. Furthermore, we denote {Z)^ := 
{x,0\WZiU\x,0). We now state the following claim. 

Claim. U computes f with probability at least p if and only if, for every x, one 
has 

(a) \{Z)x\ >2p-l and 

(b) sign{Z), = (-!)/(-). 

To prove this claim, consider the preparation of U\x,0) followed by a compu- 
tational basis measurement of the first qubit and let Px denote the probability 
that the measurement outcome is /(a;), for every a;. Then U computes / with 
probability at least p if and only if p^; > p for every x. We thus have to show 
that (a)-(b) are equivalent to the condition Px > P for all x. 

We will distinguish between the cases f{x) = and f{x) = 1. We start 
with the former case. As px and 1 — px are the probabilities of obtaining the 
measurement outcomes and 1, resp., one has 

{Z)x^Px-{l-Px)=2px-l. (11) 

Now suppose first that Px ^ P for every x (with p G (0.5, 1] as stated above). 
Conditions (a) and (b) then follow immediately. Conversely, assume that (a) 
and (b) are true. As (b) is satisfied, we have {Zi)x > and using (a) it thus 
follows that {Zi)x > 2p — 1. Invoking ([TT|) then implies that Px > P for every x, 
as desired. This completes the proof for the case f{x) = 0. 

The case f{x) = 1 is treated analogously; the main distinction is that now 
Px represents the probability of measuring 1. Consequently, dill) is replaced by 

{Z)x = (1 - Px) ~Px = l- 2px. (12) 

The remainder of the argument is completely analogous. 

5.2 Step 2 

Conditions (a)-(b) imply that the study of matchgate-computable functions re- 
duces to the investigation of matrix elements of the form (x,0|[/^ZiC/|x,0), 
where U is an arbitrary matchgate circuit. Next we investigate the most gen- 
eral form which such matrix elements may take. A complete characterization of 
this problem is obtained in the following theorem. 

Theorem 5. Let U be an n-qubit unitary operator. If U is a matchgate circuit 
then there exists an a E R" with \\a\\i < 1 such that 

{x\U^ZiU\x) =a^x (13) 

for every n-bit string x. Conversely, for every a G M" with \\a\\i < 1 there exists 
an n-qubit matchgate circuit U such that I113\) holds. 
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In the proof of this theorem we will need the following elementary fact. 

Lemma 1. Consider a vector a £ R" with \\a\\i < 1. Then there exist u, w £ R" 
with \\u\\2 = 1 and \\v\\2 < 1 such that a^ '■— u^Vk for every k G [n]. 

Proof: Define u and v by 

Uk := V|afc|/||a||i and Wfc := \/|afe| ■ ||a[|i ■ sign(afc), (14) 

for every k e [n], respectively. Obviously, Ofc = u^Vk- Moreover, ||u||2 = 1 and 
Iklb — \\a\\i < 1, as can be easily verified. This proves the lemma. D 

Proof of theorem O' We first prove the forward direction. Denote O := 
WZiU and let R be the SO{2n) rotation associated to U via (j4]). Letting p 
and p' denote the first, resp. second, row of R and using that Zi = — iciC2, it 
follows that 

O = -i[U^ciU][U^C2U] = ^p^p'.hic^c,], (15) 

where the sum is over all p,v ^ [2n]. As {x\0\x) is a diagonal entry of O 
for every x, we only need to focus on the diagonal part of this operator i.e. 
diag(O) :— '^.j.{x\0\x)\x){x\. Using the exphcit representation (jSj of the c^, it 
is easily verified that 

li fi = v 
,. ,_. ._, ^K if (^i, J^) = (2fc - l,2fc) for some fc e [n] , , 

mag^ ic^c,} -\ 7. if (^^ y) ^ (2fc, 2k - 1) for some k e [n] ^^^ 

otherwise. 

Setting Ofc := p2fc-iP2fe ~ P2kP2k~i ^r every k e [n] and using that X)^ P^pj^ is 
zero since p and p' are two distinct rows of an orthogonal matrix, it follows that 




diag(0)-^afeZfe. (17) 

k=l 

The {x, x) diagonal entry of O thus reads: 

{x\0\x)=J2'^k{x\Zk\x)^a'^x. (18) 

This shows that that (fT3|) is satisfied. Note also that |(a;|0|a::)| < 1 for every 
X since O is unitary. Moreover, it is easily verified that there exists an x such 
that a'^x = \\a\\i. This shows that ||a||i < 1. 

To prove the reverse direction, consider an arbitrary a G R" with ||a||i < 1. 
Due to lemma [1] there exist n-dimensional real vectors u and v with ||m||2 = 1 
and ||w||2 < 1 such that Ok = UkVk for every k G [n]. Furthermore, it is 
elementary that for every such u and v there exists a vector w G R" that is 
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orthogonal to u and that satisfies ||w||2 = 1 — Iklli- Consequently, the 2n- 
dimcnsional vectors 

/9 := (Ul,0,-U2,0,...,-U„,0) and p' -.^ {wi,Vi,W2,V2, . ■ . ,Wn,Vn) (19) 

are unit vectors (w.r.t. the 2-norni) that are orthogonal. Let R be any SO{2n) 
rotation having p and p' as first, resp. second row. Let U be the n-qubit match- 
gate circuit associated to R. Analogous to the proof of the forward direction of 
the theorem, a direct calculation shows that the diagonal part of C/^Zit/ equals 

n 

^[P2k-lP2k - P2kp'2k-l\Zk = X! ^kZk- (20) 

fe=l 

It immediately follows that ([T^ is satisfied. This completes the proof. D 



5.3 Step 3 

We now show how the above result leads to the proof of Theorem [TJ Invoking 
Step 1, it suffices to show that, for every n-bit boolean function / andp G (0.5, 1], 
the following arc equivalent: 

• There exists an rn-qubit matchgate circuit C/, for some m > n, such that 
conditions (a) and (b) in Step 1 hold; 

• / is a linear threshold gate with margin e > 2p — 1. 

We first prove the forward direction of the claim. Let [/ be a matchgate circuit 
on m qubits such that (a)-(b) hold. We first invoke theorem [5l This yields a 
vector a G R™ with 1-norm at most 1 such that (IT^ holds. With the notations 



a :— (ai, . . . , a„) and b :— X]i=n+i "^^ ^^ then immediately follows that 

{x,0\U''ZiU\x,0)^a'^x + b (21) 

for every n-bit string x. Together with condition (b) in Step 1 this shows that 
/ is a linear threshold gate with representation (a, &). Normalizing the vector 
(a, b) w.r.t. the 1-norm yields a normalized representation w :— 7a and 9 := 76, 
where 7"^ := ||a||i + \b\. Note that 7-1 < ||a||i < 1. Using Eq. ^ and 
condition (a) of Step 1 it follows that 

|a^x-|-6| > 2p-l. (22) 

Moreover, using the definitions of w and 6 and the fact that 7 > 1, it finally 
follows that 

\w'^x + 9\>-f{2p~l)>2p-l, (23) 

for every x. This shows that the margin of / is at least 2p — 1. 

We now show the converse. Consider an n-bit LTG /. Let e £ (0, 1] be 
its margin with associated (normalized) representation (w, 9) and denote p := 
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(e + l)/2. Due to theorem [SJ there exists an {n + l)-qubit matchgate chcuit U 
such that 

{x,y\U^ZiU\x,y) =w'^x + ey (24) 

for every n-bit string x and for every y — 0, 1. This readily imphes that 

sign((x, 0\U''ZiU\x, 0)) ^ sign(w'^x + 6) = (-!)■'''"=', (25) 

for every x. Moreover, as e is the margin of /, it foUows that 

|(x,0|C/tZiC/|x,0)| = \w^x + 0\ >2p-l (26) 

for every x.This shows that conditions (a)-(b) in Step 1 are fulfiUed. This 
completes the proof of theorem [1] 

Remark. — It follows from the above argument that an n-bit boolean function 
/ is computable by an m-qubit matchgate circuit with probability at least p, 
for some m > n, ii and only if / is computable with probability at least p by a 
matchgate circuit acting on at most n + 1 qubits. o 

6 Margins of threshold gates 

In this section we analyze the properties of margins of LTGs in more detail. In 
particular, the proofs of theorems [2] and |3] will follow from the considerations in 
this section. 

6.1 Large margins 

Next we investigate the subclass of linear threshold gates with large margin. 
We will in particular show that any LTG with margin strictly larger than 1/2 
is essentially trivial. The proof of theorem [5] will follow immediately from this 
property. 

Recall that, formally, a boolean function f on n bits is said to depend on its 
fc-th variable, if there exists an n-bit string x such that f{x) ^ f{x-,k), where 
x^k is the string obtained by flipping the fc-th bit of x. We can now state the 
following result: 

Lemma 2. Any linear threshold gate with margin e can depend on at most 
t — 1/e of its variables. 

Proof: Let / be an n-bit LTG and let e be the margin of / with associated 
normalized representation {w,9). We make the following claim. Claim: If / 
depends on its fc-th variable then \wk\ > £■ The proof of the lemma immediately 
follows from correctness of this claim, by using that |jw||i < 1 as {w,6) is a 
normalized representation. We now prove the claim. For simplicity but without 
loss of generality, we set fc := 1 and assume that / depends on its first variable. 
We prove that \wi\ > e. Denote w := {w2, ■ . . ,Wn). If / depends on its first 
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variable, there must exist an s G {±1}""-'^ such that w'^s± \wi \+9 have opposite 
signs. We consider the fohowing two possibilities: (a) w'^s + 6* > or (b) 
w'^s + 9 < 0. In case (a), it follows that w'^s + \wi\+9 > and hence we must 
have w'^s — jwil + 9 < 0. Since e is the margin of {w, 9), it follows that 

iv'^ s - \wi\ + 9 < -e. (27) 

Hence, \wi\ > e. Case (b) is treated analogously. D 

It follows that any linear threshold gate which depends on all its n variables 
can admit a margin of size at most l/n, and any LTG with constant margin 
can depend on at most a constant number of its input bits. What is more, if 
e is sufficiently close to 1 viz. e > 1/2 then / can depend on at most one of 
its variables. Indeed, it follows from lemma [2] that / can depend on at most 
t = €~^ < 2 input bits, i.e. i is (corresponding to a constant function) or 
1. Combining these considerations with theorem [TJ the forward direction of 
theorem [5] follows immediately: indeed, due to theorem [T] any boolean function 
which is matchgate-computable with probability p > 3/4 must be an LTG with 
margin e > 1/2, and hence of the form indicated. The proof of the converse 
direction of theorem [5] is elementary, as discussed in section [3] 

More general than theorem [5J we have actually showed that any function 
which can be computed by a matchgate circuit with success probability at least p 
is an LTG which depends on at most 0{p~^) of its input bits. As a consequence, 
whenever a constant success probability is considered (i.e. independent of the 
input size), only functions can be computed which depend on a constant number 
of input bits. Conversely, we will prove that any linear threshold gate depending 
on k of its input bits can be computed by a matchgate circuit with success 
probability of at most ^ + 0(fc~^). Therefore, any LTG that depends on all its 
n input bits can be computed by a matchgate circuit with probability at most 
^ + 0{n~^), which always lies at least polynomially close to 0.5. 

We conclude this section with an example. Let n be odd and let /maj denote 
the majority function on n bits i.e. /maj(a^) is iff the n-bit string x contains 
more zeros then ones. It can easily be verified that /maj is an LTG with nor- 
malized representation {'w,9) defined by Wk ■— n^^ for every k and 9 := 0. 
Moreover, the minimal value of \w'^x + 9\ over all x is easily shown to be n~^, 
showing that n~^ is a lower bound for the margin of /maj. As /maj depends 
on all its n variables, its margin is at most n~^ due to lemma [21 This shows 
that the margin of the majority function is precisely n^^ = emaj. Due to theo- 
rem [TJ there hence exists a matchgate circuit that computes this function with 
probability 

,.„.£^.l+_L^ (28) 

Moreover, any matchgate circuit computing the majority function can do so 
with probability at most pmaj . 
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6.2 Margins and integer representations 

Next we focus on the possible types of asymptotic behavior the margins of a 
family of LTGs may exhibit. We will in particular be interested in the dis- 
tinction between margins that are polynomially bounded and margins that are 
exponentially small in the input size. As discussed in section [21 the subclass 
of (families of) LTGs with polynomially bounded margin coincides with those 
functions that are matchgate-computable with poly-bounded error. In the fol- 
lowing we will in particular prove a simple characterization of the families of 
LTGs with poly-bounded margin. 

An important parameter of an LTG / in the present context will be the 
integer weight Lu{f) = uj oi f, defined to be the minimal 1-norni Huijli + \9\ of 
any possible integer representation of /. A family of LTGs {/„ : n = 1, 2, . . . }, 
where /„ acts on n bits, is said to have polynomial integer weight if the integer 
weight LOn of /„ scales polynomially with n. Interestingly, it turns out that the 
margin of an LTG and its integer weight are closely related concepts: 

Theorem 6. Consider an n-bit linear threshold gate with margin e and integer 
weight uj. Then 

i<^<2.!i±l. (29) 

Note that, in the case of large n, the second inequality in (|29l) simplifies to 
oj < 0{ — ). As an immediate corollary of theorem [SJ the following property 
follows: 

A family of linear threshold gates has poly-bounded margin if and only if 
it has polynomial integer weight. 

Combining this result with theorem [1] then proves theorem [31 

Proof of theorem [^' We first prove the first inequality. Let (v, ip) be an 
integer representation with \\v\\i + \(p\ = uj. Normalizing {v,(p) leads to the 
representation {w, 9) defined by w := uj'^v and 9 = uj~^(p. We claim that {w, 9) 
has margin at least a;~^, implying that e > uj~^ as desired. To prove the claim, 
note that the sign of v^x -\- 9 equals (— l)-'"'^) for every n-bit string x. This 
implies that v^x -t- (^ 7^ for any x. As the components of v and (f are integers, 
v^x + 'f is an integer as well for any x. Hence, the property v^x + ip ^ implies 
that \v'^x -|- (/?! > 1. It follows that \w'^x + 9\> uj^^ i.e. the margin of (w, 9) is 
at least w~^. 

Next we prove the second inequality. Let {w, 9) be a normalized representa- 
tion of / with margin e. Let d be a positive integer; for now d is arbitrary but we 
will fix a value later. Let w^fj, and 0^ represent the /x-th bit in the binary expan- 
siorO of Wk and 9, respectively. Now let Vk and (p be the rational numbers ob- 
tained by truncating the binary expansion of Wk, resp. 9, after the d-th bit (and 



^Any real number a G [1,-1] can be expanded in a unique way as a 
sign{a) JZlT^i '*m2~'^' where each a^ G {0, 1}. 
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keeping the same overall signs). That is, we define Vk — sign(w;/j) ^ Wfc^2~^ 
and (p = sign(6') Yl d^2~^ for every fc, where the sums run from 1 to d. Note 
that 

lkl|i + l¥'l< lklli + 1^1=1- (30) 

It further follows from the definitions of v and </? that \wk — Vk\ < 2^"^ and 
\ip-e\ <2-'^. This implies that 

\{v~wfx + {ip-e)\<{n + l)2-'^. (31) 

We now choose d to be the smallest positive integer such that {n + 1)2^'' is 
strictly smaller than e. Then, as \w'^x + 9\ > e, the quantity 

w'^x + e + {{v^wfx + {ip-9)} (32) 

must have the same sign as w^x + 9 for every x. Note that ([5^ coincides with 
v'^x + 9 for every x. This shows that {v, (p) is also a representation of /. But 
then also multiplying (v, ip) with 2'^ leads to a representation of /. The latter 
representation is integer, and moreover has l-norni at most 2*^, since the 1-norm 
of (v, (p) is at most 1 due to ((30)) . This shows that the integer weight of / is at 
most 2"^. Finally, for our choice of d one has 

2d<2-^^±l. (33) 

To see this, remark that d is defined to be the smallest integer strictly larger 
than log2 [^^-^j • But then d must satisfy 



d < log2 



n+1 



1, (34) 



£ 

which is equivalent to (|33| . This proves the second inequality in (|29|) . D 

7 An equivalent classical computer 

Theorem [1] connects the optimal success probability p of computing a linear 
threshold function / with any matchgate circuit with the margin e of this LTG. 
This is a somewhat peculiar connection, e.g. since the definition of the margin of 
an LTG a priori does not seem to have anything to do with probabilities. In order 
to understand this relation better, in this section we construct a very simple class 
of classical computers that are capable of (probabilistically) computing LTGs 
with precisely the same relation between p and e. 

To motivate this class, we re-iterate the following simple example which was 
discussed in section [51 Let /maj denote the n-hit majority function as before, 
and let e :— n~^ be its margin. Due to theorem [U there exists a matchgate 
circuit that computes this function with probability Pmaj defined in (j28p . and 
this is the optimal success probability which any matchgate circuit can achieve. 
Now consider the following elementary classical computation, consisting of the 
following steps: 
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• Choose an n-bit input string x. 

• Generate a random integer k between 1 and n. 

• Output the bit value Xk- 

We now claim that, for every x, the output of this computation is /maj(a^) with 
probability at least Pmaj- To see this, note that the probability px that the 
above procedure outputs is given by 

Px = \{k^[-\---^-^}\^ (35) 

n 

and the probability of outputting 1 is q^ = 1 — Px- If 2; contains more zeros 
than ones — i.e. if /majl^^) is zero — then the probability px is at least Pmajj as 
can easily be verified. Similarly, if x contains more ones than zeros, qx is at 
least Pmaj as well. Thus, for any x the output of the computation is /maj(a^) 
with probability at least Pmaj- In other words, the above classical computation 
computes the majority function with the same probability of success which can 
be achieved by the optimal matchgate circuit. Next we will show that every 
LTG / with margin e can be associated with a classical computation of the 
above nature. 

The relevant class of classical computations is defined as follows. Fix a 
probability distribution V := {T^k} on the set of integers from 1 to n+1, together 
with a string c of n + 1 bits. We define the weighted majority sampling (WMS) 
computation associated with V and c to consist of the following steps: 

• Choose an n-bit input string x. 

• Sample from the distribution 7^, yielding k <E [n + 1] with probability tt^. 

• li k < n then output the bit Zout '■— Xk (S Ck- If fc = n + 1 then output 

We say that an n-bit boolean function / is WMS computable with probability 
at least p if there exist V and c such that the above three-step procedure yields 
the output Zout — f{x) with probability at least p for every n-bit input x. We 
prove that the classes of WMS computable and matchgate-computable functions 
precisely coincide. 

Claim. Let p S (0.5, 1]. A function is m,atchgate-com,putahle with probability at 
least p iff this function is WMS computable with probability at least p. 

Consider a WMS computation with associated V and c. It will be convenient to 
consider slightly modified computation where the output is Zout = (— 1)^°"* = 
±1 instead of the bit Zout- This will facilitate notation in the proof (but does not 
play any essential role otherwise). Further, we let (zout)a; denote the expected 
value of Zout given that x is the input of the computation. 

We now prove the claim. Let / be an arbitrary n-bit boolean function. It 
follows from section [Ql that / is machgate-computable with probability at least 
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p iff there exists a matchgate circuit U acting on to > n qubits such that, for 
every x, one has 

(a) \{Z)x\ >2p~l and 

(b) sign(Z), = (-!)/(-). 

Furthermore, using an argument analogous to in section 15.11 it can easily be 
shown that / is WMS computable with probability at least p if and only if 
there exist V and c such that, for every x: 

(a') \{zout)x\ > 2p- 1 and 

(b') sign(zout). = (-l)-^("\ 

for every x. We thus have to prove that, for every function, conditions (a)-(b) 
hold for some matchgate circuit U iff (a')-(b') hold for some V and c. 

Suppose first that (a')-(b') are satisfied for some V and c. We define w £ K" 
and 6* G R by 

Wk := {-ly-TTk and 6 := (-l)^"+i7r„+i, (36) 

for every fee [n]. Note that the vector (zw, 9) has unit 1-norm. Now let x be an 
arbitrary bit string and run the WMS computation as described above. Then 
the expected value of iout is 

(zont). = |X^^fc(-l)"'=+''' I +^„+i(-l)^'" =u;^£ + 0. (37) 

Due to theorem [SJ there exists a matchgate circuit U on n + 1 qubits such that 

(Zi), := {x, 0\U^ZiU\x, 0)=w'^x + e^ (zout).- (38) 

As {Zi)x — {zout)x for every a;, it follows that conditions (a)-(b) are satisfied 
fort/. 

To prove the converse, consider an ?7i-qubit matchgate circuit U such that 
(a)-(b) hold. Due to theorem[5l there exists {v,(p) with 1-norm at most 1 such 
that 

{Zi):, = (a-,0|[/^ZiC/|x,0) ^v'^x + ip (39) 

for every n-bit string x. Normalizing (v, ip) w.r.t. the 1-norm yields a normalized 
representations :—^jvan(i9 ;= 79?, where 7""'^ := ||t'||i-|-|(/9|. Note that 7^"'^ < 1. 
Now choose V and c such that ([55]) is satisfied. Using an argument similar to the 
first part of the proof, the expected value of the associated WMS computation 
is {zoMt)x = w^x + 9. We thus have 

(^out)x = ^{v^x + (p) = l{Zi)^ (40) 

for every x, where 7 > 1. Using the identity (iout)x — 7(^i)x and the fact that 
(a)-(b) hold, it immediately follows that (a') and (b') are satisfied for {V,c). 
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